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00 , 

Q\ \ Laplacian Abelian Projection is discussed. This term refers to the use of a new ("Lapla- 

. cian") gauge fixing prescription for implementing the Abelian Projection of QCD. The gauge 

condition is based on the lowest-lying eigenvector of the covariant Laplacian operator in the 
, | ' adjoint representation. 

This Laplacian gauge fixing procedure is free of the ambiguities which plague lattice 
simulations which work with the popular Maximally Abelian Gauge. Furthermore, Laplacian 
gauge fixed configurations enjoy a natural kind of smoothness. These two properties are 
' crucial for a reliable determination of physical quantities using the Abelian Projection. 

We also examine a new, Higgs-field-like observable which emerges as a by-product of the 
method. This quantity can be used to identify magnetic monopoles in a way independent 
of the traditional prescription. It is argued that physically relevant magnetic monopoles 
are accomodated well by the Laplacian method, while they are suppressed (too) strongly in 
f"^ , Maximally Abelian Gauge. 

Finally, first evidence of abelian dominance in the Laplacian Abelian Projection is pre- 
£f) ' sented. 

o ' 

. §1. Panorama 

't Hooft's Abelian Projection paper & has been the point of departure for many 



simulations of lattice gauge theories in the past decade. His idea was that (partial) 
gauge fixing could expose the degrees of freedom relevant in the long-distance regime 
of QCD. He showed that if an abelian subgroup U(l) Ar_1 of the non-abelian gauge 
group SU(N) is left unfixed, one obtains an abelian gauge theory whose gauge field 
content consists of "electrically" charged vector fields and "magnetic" monopoles in 
addition to the N — 1 abelian "photons" . The magnetic monopoles arise as defects 



in the gauge fixing. 

An example of su ch a partial gauge condition in the continuum is given by the 
covariant gauge f 



E 



(d,^iAl)A± = 0. (1-1) 



In this gauge, the non-abelian components A 1 ' 2 are covariantly constant with respect 
to the abelian subgroup which is taken in the 3-direction. There is no gauge condition 
on the abelian component A 3 . 

On the lattice, the most popular gauge is the "Maximally Abelian Gauge" 
(MAG)!), which is defined as the configuration {llj^J} for which the gauge trans- 



* 5 ETH preprint SCSC-TR-98-01. Invited talk at YKIS'97 (Kyoto, Japan, December 1997). 
**•* E-mail address: arjan@scsc.ethz.ch 

***) We shall restrict ourselves to SU(2) throughout this presentation. 
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formation {fi x } minimizes the functional 



;i- 2 ) 



X,fl 



The lattice gauge field U^ x transforms as U^J 



These gauge conditions (1-1, 1-2]) are in fact related a). In the naive continuum 
limit, Eq. (|1-2|) becomes 



(4 fi)1 ) 2 + (4 ft)2 ) : 



;i-s) 



and the stationarity condition which its minimum has to satisfy is precisely Eq. (|1-1| ) . 

Many interesting results have been obtained using the MAG, as we will briefly 
review in Sect. |2[ This cannot be said of most other Abelian Projection gauges. The 
reason is presumably that the MAG favours relatively smooth gauge fields, unlike 
these other gauges. There are a number of important drawbacks associated to the 
MAG, however. 

Most importantly, it is impossible to implement the gauge condition unambigu- 
ously, in practice. This means, for example, that the result depends on the gauge 
in which one happens to start the gauge fixing algorithm. In other words, gauge 
invar iance is lost. 

Secondly, magnetic monopoles are not identified according to 't Hooft's prescrip- 
tion, by looking for singularities in the gauge fixing (in fact one never attempts to 
find those singularities in the MAG!), but instead using the U(l) procedure based 
on elementary cubes a. It is not clear if there exists an exact relation between the 
monopoles in the two definitions. 

The Laplacian method discussed here manifestly solves the first problem. In 
addition, it seems to induce a natural solution to the second drawback: we find 
evidence that the gauge fixing singularities are correlated with the "elementary-cube 
monopoles" in this gauge, and that they signal the presence of physically relevant 
monopoles. 



§2. Silver Pavillion: Maximally Abelian Gauge 



Let us briefly sketch what lattice simulations using the MAG have achieved in 
recent years. 

What one does in lattice simulations, is work with the abelian fields after the 
gauge fixing only. So one generates configurations according to the full SU(2) Haar 
measure, then one puts these configurations in the chosen gauge, and subsequently 
one throws away the non-abelian components A 1 ' 2 : all observables are constructed 
using the abelian field A 3 only. 

For example, one can construct "abelian Wilson loops" in this way. Surprisingly, 
one finds evidence W for the phenomenon of "abelian dominance" hypothesized in 
Ref. 6: these abelian Wilson loops show confining behaviour with a string tension 
compatible with the one obtained from the full SU(2) gauge fields. 
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Later it was found furthermore that Wilson loops made of the magnetic monopole 
fields contained in this abelian component reproduce the string tension as well 
("monopole dominance" ) Q) . Issues related to chiral symmetry breaking^ have also 
been studied using the abelian fields and the monopoies in the MAGeP. Let us 
emphasize that these monopoies are defined using the elementary-cube definition of 
Ref. 4, not by identifying singularities in the gauge fixing procedure. 

More recently, studies of monopole effective actions and their renormalization 
group behaviour have yielded very interesting results™. 

Other gauge conditions for the Abelian Projection have not shown convincing 
evidence of abelian or monopole dominance (and monopole effective action results 
are not available for other gauges than MAG) . The reason is presumably that MAG is 
a smooth gauge, while the others are not. By "smooth" we mean, roughly speaking, 
that (the non-abelian parts of) the link matrices are close to unity, and do not vary 
too wildly from one site to the next. This kind of smoothness is required for a reliable 
mapping {U^ :X } — > {A^ jX } of lattice to continuum gauge fields, which is inherent in 
the projection procedure implemented in lattice simulations. 

The successes of the MAG are overshadowed by an important shortcoming, how- 
ever: It is ambiguous! The gauge condition is implemented by means of a local itera- 
tive minimization algorithm, which often ends up in one of a number of local minima 
of the functional ([L^) • This means that two gauge equivalent starting configurations 
may lead to different gauge fixed configurations. In other words: gauge invariance is 
lost, and the physics extracted from the abelian projection using the MAG becomes 
unreliable. These local minima are sometimes called "Gribov copies". However, 
most of them are "lattice Gribov copies", UV artefacts of the lattice procedure, 
which should be distinguished from "true" Gribov copies. 

Another weak point is that the magnetic monopoies are not detected by looking 
for singularities in the gauge fixing procedure in the spirit of 't Hooftu), but instead 
by projecting out the abelian component and constructing monopoies according to 
the elementary-cube construction l§ familiar from U(l) theories. It is not clear that 
these are always related. One can consider the physical example of a dyon, to 
gain some insight into this relation (see Ref. 10 for a discussion of the dyon in 
different abelian projection gauges). Here one finds that both definitions point 
towards a monopole singularity at the origin. The singularities of the MAG fixing 
are found to occur when the non-abelian fields diverge (become of order 0(a^ 1 ) 
on the lattice), while the cube definition looks for Dirac monopoies: points where 
the abelian component blows up. For the dyon, which is a regular solution, these 
singularities in the abelian and non-abelian gauge field components have to coincide, 
but for a general elementary-cube Dirac monopole found in simulations it is not clear 
that it comes from a physically relevant monopole. 

In fact, one can argue that physical monopoies are suppressed too strongly by 
the MAG. For the (continuum) dyon in MAG (IT), we know the non-abelian fields 
A 1,2 blow up, so the functional ( |T3j ) must accomodate a large contribution from the 
central region of the dyon. This goes against its tendency to minimize this quantity. 
So here we have another drawback of the MAG: it unnaturally attempts to suppress 
physical monopoies. 
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The Laplacian method for Abelian Projection will be seen to kill all these birds 
with one stone. It is unambiguous, yields smooth gauge fixed configurations, does not 
unduly suppress physical monopoles, and there is evidence that the elementary-cube 
monopoles have physical significance. 

2.1. MAG as a spin (glass) model 

To set the scene for the presentation of the Laplacian Abelian Projection in 
Sect. H|, let us try to get some insight into what the MAG actually doesll2P 
Going back to the functional S of Eq. (|1-2|) , we note that it can be written as 



SUM) = E I 1 - ^ "'■rl' i ,r'l>r. il l ■„.,■} (24) 

= E{i-£#*?MUa}. (2-2) 

x,n a ,b 

where 

R% = ±Tr [a a U^ x a b U+ x ] (2-3) 

is the lattice gauge field matrix in the adjoint representation. Here <P is defined as 

3 

$ x = f2+a 3 n x = Efc- (24) 

a=l 

It satisfies the constraint 

E(^) 2 = l- (2-5) 

a=l 

4> is a very natural variable to parametrize the gauge transformation Q: the 
original gauge transformation matrix Q x G SU(2) is determined up to the residual 
U(l) gauge invariance only, but <1> X E S 2 ~ SU(2)/U(1) is invariant under f2 x — > 
V X Q X for diagonal matrices V x G 17(1) C SU(2): 

$ x = f2+a 3 n x Q+V+a 3 V x n x = $ x . (2-6) 

(In other words, <P is neutral with respect to the residual U(l) symmetry.) 
Eq. (2-2) is the lattice discretization of the continuum functional 

S A (<f>) = [ \{D^)\ (2-7) 



,v 2 

which could also have been obtained directly from Eq. ( |l-3p . This is the covariant 
kinetic action of a unit-length adjoint spin field in the background of the gauge field 
configuration which one wants to gauge fix. It is this functional which the MAG 
attempts to minimize. 

Because of the quasi-random nearest-neighbour couplings between the spins in 
this action (2-2), it is similar to a spin glass model. It is the length-1 constraint on 
the spin variables which is responsible for the difficulties in solving this minimization 
problem. 
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§3. Golden Pavillion: Laplacian Abelian Projection 

The Laplacian method is based on relaxing the length- 1 constraint on the 
<j> fields in the minimization problem (2-2 j2-7|) of the MAG. 

In this case the "spin glass action" (2-2 j2-7| ) becomes an unconstrained quadratic 
form which is minimized by the the lowest-lying eigenvector of the covariant Lapla- 
cian operator — = — ^ D^(i?)Z)^(J2). (We do of course maintain an overall 
normalization on the eigenvectors.) On the lattice, —□(/?) is the 3V x 3V real sym- 
metric matrix 

Its eigenvalues and eigenvectors can be computed numerically in an unambiguous way 
and to arbitrary precision. One can use the Lanczos method (which in addition to 
eigenvalues can provide the corresponding eigenvectors), or a Rayleigh-Ritz method, 
for example. 

The only situation in which an ambiguity might arise is when the lowest eigen- 
value is degenerate. The two (or more) corresponding eigenvectors would lead to 
different (but all smooth) gauge field configurations, both satisfying the gauge con- 
dition exactly. One would have a true, physical Gribov copy (as opposed to the 
"lattice Gribov copies" afflicting the MAG). However, this set of true Gribov copies 
has measure zero, and they are never encountered in practice: The two lowest eigen- 
values are always found to be separated by an amount which is orders of magnitude 
larger than the numerical precision with which they can be resolved. 

It remains to be discussed how the lowest mode of the Laplacian determines the 
gauge fixing. By relaxing the length- 1 constraints on the fields <j> x we have of course 
abandoned their straightforward interpretation in terms of gauge transformation 



matrices £2 X through Eq. (2-4). 

There are two equivalent ways to see how the lowest eigenvector 0" defines the 
required gauge transformation Q x . One way is to normalize 4> x site by site, by writing 



PxFx , (3-2) 




Pz=(2>")l > ( 3 - 3 ) 

Va=l 

and to use the normalized vector to obtain Q x (up to its U(l) arbitrariness of course): 

3 

f2+* 3 f2 x = ^fr. (3-4) 

a=l 

The other way is to realize that the matrix which rotates the 3-vector <j) x into the 
abelian direction (which we chose to be the 3-direction), is nothing but the required 
gauge transformation fl x in the adjoint representation. In this picture, the U(l) 
arbitrariness corresponds to rotations about the direction of <p x . 

*•* This follows a similar idea in the context of Landau gauge fixing on the latticeB. 
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This procedure is ill-defined only when p x = for some x. At such points there 
is no way to decide which Q x to take. In other words, since \(f)x\ — 0? the symmetry is 
not broken to U(l) but one retains the full SU(2). This is a singularity in the gauge 
fixing procedure, of the type that identifies a magnetic monopoleEP. Note that this is 
very similar to the zero of the Higgs field in the true 't Hooft-Polyakov monopole in 
the Georgi-Glashow model! This is a promising feature, suggesting a direct relation 
of these gauge fixing singularities with physical monopoles. We shall discuss this in 
detail in Sect. 3.2. 

3.1. Laplacian Abelian dominance 

We will present preliminary results concerning abelian dominance using the 
Laplacian Abelian Projection. 

Fig. 1 shows the "square" Creutz ratios 



X (L + l/2,L + l/2) 



In 



(W(L,L)){W(L + 1,L + 1)) 
(W(L,L + 1))(W(L + 1,L)) 



(3-5) 



(W(R,T) denotes an R x T Wilson loop) plotted against L + 1/2, for two SU(2) 
simulations. Shown are the full SU(2) Creutz ratio and the "abelian Creutz ratios" , 
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Fig. 1. Abelian dominance: comparison of full and abelian "square" Creutz ratios for two SU(2) 
simulations. The Laplacian method (LAG) seems to be slightly "more abelian dominant" than 
the Maximally Abelian Gauge (MAG). 
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formed from the abelian gauge field components after abelian projection using the 
MAG and the Laplacian Abelian Gauge (LAG). The errors (not shown) on the data 
points at L + 1/2 = 4.5 are large, especially for the 16 4 simulation at (3 = 2.5. 

One observes that xlag is about the same as xmag- Both xlag an d Xmag 
approach XFuii at large distance, but they are smaller than XFuii at short distances. 
This is in agreement with the fact that the smaller Coulomb term in the abelian 
potential leads to a smaller slope dV/dr in that regime. Thus we find that the 
LAG shows abelian dominance just like the MAG. In fact, the data seem to suggest 
that xlag m ay even lie slightly closer to XPuii than xmag- However, more accurate 
calculations are needed before this can be concluded with confidence. 

The 12 4 graph shows in addition the abelian Creutz ratios for gauge fixing to 
the MAG after preconditioning the configurations with the LAG. They seem to 
be roughly equal to xmag as calculated without LAG-preconditioning, and again 
slightly smaller than xlag- (This gives some indication of the size of ambiguities in 
MAG gauge fixing.) 

Hence we find first evidence that the LAG method passes the test of abelian 
dominance; more detailed computations as well as a study of monopole dominance 
are under way. 

3.2. The Higgs-like degree of freedom p 

We want to examine in some more detail the origin and physical relevance of the 
Higgs-like field and in particular its length p x , which emerges in the Laplacian 
gauge fixing procedure. 

First of all, note that it is a gauge-invariant property of the configuration, in the 
sense that (unlike the lowest eigenvector (f> x itself) it does not depend on the initial 
gauge of the full SU(2) configuration. 

We found that zeroes in the field <fi correspond to singular points in the gauge 
fixing, and hence to monopoles in 't Hooft's Abelian Projection picture. This is remi- 
niscent of the Georgi-Glashow model, where 3-dimensional vortex-like configurations 
of the Higgs field (with zeroes in their centre) are characteristic of the presence of 
a 't Hooft-Polyakov (tHP) magnetic monopole. One wonders how far this analogy 
can be extended to the Higgs-like field that emerges in the Laplacian gauge (LAG). 

Let us consider what form a classical dyon takes in the LAG. The answer can 
be obtained by the following argument. The LAG attempts to minimize Sa(4>) = 
/ (D^cp) 2 in the background of the gauge field configuration of the classical dyon. 
The spatial gauge field components of this dyon equal those of the tHP monopole 
in the Georgi-Glashow model, while their th or 4 th components differ. If we forget 
this latter difference for a moment, we conclude that since Sa(4>) has exactly the 
form of the Higgs kinetic action in the Georgi-Glashow model, the static Higgs field 
part of the tHP monopole solution in the BPS limit is exactly what the Laplacian 
eigenvector (j) x should look like for the dyon in the LAG. Finally, it is clear that 
the non-zero 4 th component of the dyon solution does not change this conclusion: it 
(also) has the form of the static Higgs field of the tHP monopole in the BPS limit, 
so its contribution to (D^cj)) 2 is zero. 

So the lowest eigenvector of the covariant Laplacian in the background of a 
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classical dyon of scale p, in radial gauge is given by 

4>% = 5 ak x k p{r) , (3-6) 
cosher 1 

P( r ) = M^-r • (3-7) 

smn pr r 

This implies first of all that the Laplacian gauge is reached by rotating the Higgs 
field into the 3-direction. This is simply unitary gauge, and the 3 rd component of 
the gauge field configuration will have the form of a Dirac monopole field. (So the 
elementary-cube procedure for detecting monopoles will identify it, just like it does 
in the MAG.) At the same time, we see that, in Laplacian gauge, the classical dyon 
shows up as a region (of extent p^ 1 ) of small p, with a zero in its centre. 

This also allows us to deduce that the dyon is not suppressed in an unnaturally 
strong fashion in the LAG, as opposed to the MAG case (cf. the discussion in Sect. 
||). The point is that the quantity that is minimized in the LAG is not Sjj{Q) of Eq. 

€3) but 

Su{(2,p) = 2^> x . [(u 1 ^) 2 + (n^) 2 )] p x+fl (3-8) 

(with an overall normalization imposed on p). Here we have written the link matrices 
U^ x as u x ^ + iu x ^a a . In the continuum formulation, this corresponds to 

' ' 2 



S A {Q,p) = -^p 2 



(3-9) 



The value of p tends to zero as the non-abelian field components A 1 ' 2 blow up when 
the singularity is approached, thus moderating the contribution of the dyon to the 
functional fl3-9|). So the LAG automatically arranges things such that dyons are 
accomodated in the Laplacian gauge without a large penalty. 

As a matter of fact, this may also explain the somewhat higher monopole den- 
sity found using the LAGO: physical monopoles (dyons) are not too vehemently 
suppressed, unlike what is the case for the MAG. 

This leads us to the following question: what is the physical significance of the 
magnetic monopoles found using the conventional elementary-cube procedureEr, in 
LAG and MAG? 

For the LAG, it is natural to ask, in view of the above, whether regions of small p 
(implying the presence of physical monopoles, identified as gauge fixing singularities 
in the Laplacian Abelian Projection) are correlated with the location of monopoles 
defined by means of the elementary-cube procedure. This correlation was visualized 
in the video shown at the workshop, and included in the CD-ROM version of these 
Proceedingsliij' . For the MAG, we studied a similar correlation. 

Fig. 2 shows one frame from this video. The two panels correspond to one and 
the same situation, visualized in different ways to highlight different features. 

The main observations are: 

i) The regions of small p (smaller than a cutoff of 0.4, relative to an average of 
1.0) form large-scale structures, rather than looking like widely scattered UV noise. 

ii) There is a clear correlation between these regions and the location of mono- 
poles identified using the elementary-cube method. In Fig. 2 this is particularly clear 
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Fig. 2. Correlation between small-p regions and elementary-cube monopoles for Laplacian Abelian 
Projection. The two panels are different visualizations of the same situation, taken from a 12 a x 4 
simulation at /3 = 3.0, with periodic boundary conditions. The box shows the 12 x 12 x 4 x-y-t 
hyperplane at some (dual) value of z. The blue regions correspond to small p (in the bottom 
panel, only some of the boundaries are shown). The orange tubes represent magnetic monopole 
currents running in the x-y-t hyperplane, the small yellow and red balls are monopole currents 
in the 2-direction orthogonal to this hyperplane. 
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in the bottom panel, where one sees that the monopoles are enclosed in the small-p 
regions. 

iii) Furthermore, although somewhat unrelated to the present discussion, it is 
observed that (in the high-temperature phase) both small-p regions and monopole 
loops appear to align along the euclidean time direction. This is in agreement with 
the prediction of Ref. 3, which describes the QCD vacuum as a condensate of dyon- 
like configurations. Should this feature survive in full QCD (cf. the conclusions 
of Ref. JJ5), then it would be at odds with expectations based on instanton liquid 
models Ej' which predict instanton-antiinstanton molecules (dipoles) aligned along 
the euclidean time direction in the full theory. 

iv) For comparison, we looked for a similar correlation in the case of the MAG. 
Since we do not have a quantity like p at our disposal in this gauge, we considered 
instead the correlation of the monopole locations with the "non-abelian residue", 
i.e., the local value of the integrand of ( |l-2j ). The behaviour is strikingly different. 
In contrast with what is observed in the LAG, this non-abelian residue looks like 
purely short-range UV noise (countless small blue "splinters" as opposed to the few 
large blue blobs in Fig. 2), without any long-distance structure. No correlation with 
the monopole locations is observed either. 

These observations lead to the conclusion that the p field signals important 
physical structure, and that the monopoles found in the LAG have physical meaning. 
No such evidence is available for the MAG. 

§4. Meditation 

We have presented a new, "Laplacian" method (LAG) to perform the Abelian 
Projection on the lattice, which has several advantages over the popular Maximally 
Abelian Gauge (MAG). 

Most importantly, it is not afflicted with the ambiguities inherent to the local 
iterative minimization procedure needed in the MAG. It is a manifestly unambigu- 
ous procedure, which leads to gauge fixed configurations endowed with a natural 
(Laplacian) kind of smoothness. 

In addition, it treats physical monopoles mildly, contrary to the MAG. The large 
penalty which they receive in the MAG and which gives rise to a strong suppression, 
is compensated by the Higgs-like field in the LAG, which automatically arranges 
itself so as to accomodate the monopoles comfortably. We conclude that the MAG 
suppresses physical monopoles too strongly. Or in other words, it ignores important 
physical information which is contained in the variable p. 

Finally, there is evidence that, in the LAG, the monopoles found using the 
elementary-cube method are correlated to those defined by the singularities in the 
partial gauge fixing. This evidence was obtained from an analysis of the classical 
dyon on the one hand, and from studying the correlations between the Higgs-like 
variable p and magnetic monopoles in simulations on the other hand. 

Furthermore, we found first evidence for abelian dominance using the Laplacian 
Abelian Projection, similar to the MAG. This study will be pursued in more depth 
and extended to monopole dominance in the near future. 
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